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Abstract 

The process e"'"e~ — >• Z77 — > qq77 is studied in 0.5 fb~^ of data collected with 
the L3 detector at centre-of-mass energies between 130.1 GeV and 201.7 GeV. Cross 
sections are measured and found to be consistent with the Standard Model expec- 
tations. The study of the least energetic photon constrains the quartic gauge boson 
couphngs to -0.008 GeV'^ < ao/A^ < 0.005 GeV-^ and -0.007 GeV"^ < aJA^ < 
0.011 GeV-^ at 95% confidence level. 



Submitted to Phys. Lett. B 



1 Introduction 



The LEP data offer new insight into the Standard Model of electroweak interactions fl] by 
investigating the production of three gauge bosons. Results were recently reported on studies 
of the reactions e+e~ — >• Z77 |^ and e+e~ W"'"W~7 [Sllll- This letter describes the extension 
of the study of the e"'"e^ Z77 process to centre-of-mass energies, a/s, between 130 and 202 
GeV. Final states with hadrons and isolated photons are considered to select Z77 qq77 
events. 

In the Standard Model, the e"^e~ — > Z77 process occurs via radiation of photons from the 
incoming electron and/or positron. One possible diagram is presented in Figure^. 

The e"'"e^ Z77 signal is defined by phase-space requirements on the energies and 
angles 9^ of the two photons, and on the propagator mass V^: 

E^>5 GeV (1) 

|cos0^|<O.97 (2) 
\V?-mz\<2Tz (3) 

where mz and Fz are the Z boson mass and width. In the following, hadronic decays of the Z 
boson are considered. Events with hadrons and initial state photons falling outside the signal 
definition cuts are referred to as "non-resonant" background. 

A single initial state radiation photon can also lower the effective centre-of-mass energy 
of the e"'"e~ collision to mz, with the subsequent production of a quark- ant iquark pair. This 
photon can be mistaken for the most energetic photon of the e"'"e~ — > Z77 —>■ qq77 process. Two 
sources can then mimic the least energetic photon: either the direct radiation of photons from 
the quarks or photons originating from hadronic decays, misidentified electrons or unresolved 
7r°s. These background processes are depicted in Figures [TJo andHt, respectively. 

In order to compare experimental results with e"'"e~ qq77 matrix element calculations, 
a further requirement is applied on the angle between the photons and the nearest quark: 

cos ^T,q < 0.98. (4) 

This cut avoids coUinear divergences. Its inclusion makes the signal definition used here different 
from the previous one Signal cross sections calculated with the KK2f Monte Carlo program 
jS] range from 0.9 pb at = 130.1 GeV down to 0.3 pb at = 201.7 GeV. 

The Z77 final state could also originate from the s-channel exchange of a Z boson, as 
presented in Figure [Hi. This process is forbidden at tree level in the Standard Model, but 
it is expected to occur in the presence of Quartic Gauge boson Couplings (QGC) beyond the 
Standard Model. 

2 Data and Monte Carlo Samples 

This measurement uses data collected with the L3 detector at LEP in the years from 1995 
through 1999, at centre-of-mass energies between y/s = 130.1 GeV and a/s = 201.7 GeV, for 
a total integrated luminosity of 0.5 fb~^. The centre-of-mass energies and the corresponding 
integrated luminosities are listed in Table Q Given their relatively low luminosities, the y/s = 
130.1 GeV and ^/s = 136.1 GeV data sample are combined into a single luminosity averaged 
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sample at y/s = 133.1 GeV. Similarly the -y/i = 161.3 GeV and y/s = 172.3 GeV samples are 
merged into a single sample at ^/s = 166.8 GeV. 

The KK2f Monte Carlo program is used to generate e"'"e~ qq(77) events, that are assigned 
to the signal or the background according to the criteria (I) — ©. The hadronisation process 
is simulated with the JETSET ^ program. Other background processes are generated with 
the Monte Carlo programs PYTHIA [Zj (e+e" ^ Ze+e" and e+e" ^ ZZ), KORALZ ^ 
(e+e- ^ r+r-(7)), PHOJET P (e+e" ^ e+e" hadrons) and KORALW Hm for W+W" 
production except for the ez/eqq' final states, generated with EXCALIBUR [TT] . 

The L3 detector response is simulated using the GEANT [12j and GHEISHA ^S] programs, 
which model the effects of energy loss, multiple scattering and showering in the detector. Time 
dependent detector inefficiencies, as monitored during data taking periods, are also simulated. 





Integrated Luminosity 


(GeV) 


(pb-i) 


133.1 


12.0 


166.8 


21.1 


182.7 


55.3 


188.7 


176.3 


191.6 


29.4 


195.5 


83.7 


199.5 


82.8 


201.7 


37.0 



Table 1: Average centre-of-mass energies and corresponding integrated luminosities of the data 
samples used for this analysis. 



3 Event Selection 

The e"'"e~ — Z77 qq77 selection demands balanced hadronic events with two isolated 
photons and small energy deposition at low polar angle. Selection criteria on photon ener- 
gies and angles follow directly from the signal definition as -E^ > 5 GeV and | cos6'-y| < 0.97. 
The invariant mass Mqq of the reconstructed hadronic system, forced into two jets using the 
DURHAM algorithm ^1], is required to be consistent with a Z boson decaying into hadrons, 
72 GeV < Mqq < 116 GeV. 

The main background after these requirements is due to the "non-resonant" production 
of two photons and a hadronic system. The relativistic velocity jSz = Pz/Ez of the system 
recoiling against the photons, calculated assuming its mass to be the nominal Z mass, is larger 
for part of these background events than for the signal and an upper cut is used to reject those 
events. It is optimised for each centre-of-mass energy, as listed in Table El 

Other classes of background events, shown in Figure lb and Figure Ic, are rejected by an 
upper bound on the energy E^i of the most energetic photon. This requirement, presented in 
Table suppresses the resonant return to the Z, whose photons are harder than the signal ones. 
A lower bound of 17° on the angle u between the least energetic photon and the closest jet is also 
imposed. This requirement is more restrictive than the similar cut on cos6'^q included in the 
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signal definition. Data and Monte Carlo distributions of the selection variables are presented 
in Figure El for the data collected at ^/s = 192 GeV— 202 GeV when selection criteria on all the 
other variables are applied. Good agreement between data and Monte Carlo is observed. 



x/i (GeV) 


133.1 


166.8 


182.7 


188.7 


191.6 


195.5 


199.6 


201.7 


/3z< 
E^^{ GeV) < 


0.48 
31.9 


0.61 
55.0 


0.64 
67.6 


0.66 
69.8 


0.66 
72.8 


0.67 
74.2 


0.69 
75.8 


0.70 
76.6 



Table 2: Energy dependent criteria for the selection of e"'"e Z77 qq77 events. 



4 Results 

The signal efficiencies and the numbers of events selected in the data and Monte Carlo samples 
are summarised in Table 3. The dominant background is hadronic events with photons. About 
half of these are "non-resonant" events. In the remaining cases, they originate either from final 
state radiation or are fake photons. 

A clear signal structure is observed in the spectra of the recoil mass to the two photons, 
as presented in Figure El for the = 192 GeV — 202 GeV data sample and for the total one. 
The e"'"e~ — *■ Z77 —>■ qq77 cross sections, a, are determined from a fit to the corresponding 
spectra at each ^/s. Background predictions are fixed in the fit. The results are listed in Table 
13] with their statistical and systematic uncertainties. The systematic uncertainties on the cross 
section measurement are of the order of 10% [2]. The main contributions arise from the signal 
and background Monte Carlo statistics (6%) and a variation of ±2% of the energy scale of the 
hadronic calorimeter (6%). A variation of ±0.5% of the energy scale of the electromagnetic 
calorimeter does not yield sizable effects. Other sources of systematic uncertainties are the 
selection procedure (3%) and the background normalisation (3%). The latter is estimated 
by varying by 10% the normalisation of the "non-resonant" background, as estimated from a 
comparison between the KK2f and PYTHIA Monte Carlo predictions for hadronic events with 
photons, and by 20% that of the other backgrounds. Uncertainties on the determination of the 
integrated luminosity are negligible. 

The measurements are in good agreement with the theoretical predictions a^^, as calculated 
with the KK2f Monte Carlo program, listed in Table El The error on the predictions (1.5%) 
is the quadratic sum of the theory uncertainty j^] and the statistical uncertainty of the Monte 
Carlo sample generated for the calculation. These results are presented in Figure |31 together 
with the expected evolution with y/s of the Standard Model cross section. 

The distribution of the recoil mass to the two photons for the full data sample, presented 
in Figure Eb, is fitted to calculate the ratio -RZ77 between all the observed data and the signal 
expectation. The background predictions are fixed in the fit, which yields: 

i?Z77 = ^ = 0.85 ±0.11 ±0.06 

in agreement with the Standard Model. The first uncertainty is statistical while the second 
is systematic. The correlation of the energy scale and background normalisation uncertainties 
between data samples is taken into account. 
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1^0 + 08 


11 5 + 06 

± ± . _1_ W . VJ 


3 4 + 05 

O . t: _1_ W . 


1 3 + 05 


201.7 


48 


9 


6.9 ± 0.5 


5.2 ± 0.4 


1.7 ± 0.3 


0.06 ± 0.02 



Table 3: Yields of the e"'"e~ — > Z77 — > qq77 selection. The signal efficiencies e are 
given, together with the observed and expected numbers of events. The right half 
of the table details the composition of the Monte Carlo samples with A^^ denoting 
the signal, A*"^^ the qq and N^^^^"^ the other backgrounds. The uncertainties are 
statistical only. 



(GeV) 


a (pb) 


a^^ (pb) 


133.1 


0.70 + 0.40 + 0.07 


0.923 + 0.012 


166.8 


0.17 + 0.13 + 0.02 


0.475 + 0.006 


182.7 


0.36 + 0.13 + 0.04 


0.379 + 0.004 


188.7 


0.34 + 0.06 + 0.03 


0.350 + 0.004 


191.6 


0.09 + 0.09 + 0.01 


0.326 + 0.004 


195.5 


0.30 + 0.11 + 0.03 


0.321 + 0.004 


199.6 


0.28 + 0.11 + 0.03 


0.304 + 0.004 


201.7 


0.50 + 0.18 + 0.05 


0.296 + 0.003 



Table 4: Results of the measurements of the e+e~ Z77 — > qq77 cross section, a, 
with statistical and systematic uncertainties. The predicted values of cross sections, 
cr^^, are also listed. 



5 Study of Quartic Gauge Boson Couplings 

The contribution of anomalous QGCs to Z77 production is described by two additional dimension- 
six terms in the electroweak Lagrangian jTHtlTH]: 

TlOl -> -> 

where a is the fine structure constant, F^^ is the field strength tensor of the photon and 
is the weak boson field. The parameters oq and Oc describe the strength of the QGCs and A 
represents the unknown scale of the New Physics responsible for the anomalous contributions. 
In the Standard Model, = Oc = 0. A more detailed description of QGCs has recently 
appeared [T7]. Indirect limits on QGCs were derived from precision measurements at the Z 
pole [IH!. 

Anomalous values of QGCs are expected to manifest themselves via deviations in the total 
e"'"e~ — > Z77 cross section, as presented in Figured In the Standard Model, Z77 production 
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occurs via bremsstrahlung with the low energy photon preferentially produced close to the 
beam direction. The QGC s-channel production results instead in a harder energy spectrum 
and a more central angular distribution of the least energetic photon [T^. Distributions for 
this photon of the reconstructed energy, the cosine of the polar angle and the transverse mo- 
mentum for the full data sample are compared in Figure El with the predictions from signal 
and background Monte Carlo. Predictions in the case of a non zero value of oq/A^ or ac/A^ 
are also shown. They are obtained by reweighting j2] the Standard Model signal Monte Carlo 
events with an analytical calculation of the QGC matrix element Monte Carlo studies 

indicate the transverse momentum as the most sensitive distribution to possible anomalous 
QGC contributions. A fit to this distribution is performed for each data sample, leaving one 
of the two QGCs free at a time and fixing the other to zero. It yields the 68% confidence level 
results: 

ao/A^ = -0.002^°:°°^ GeV-2 and aJK^ = -0.0011°;°°^ GeV'^ , 

in agreement with the expected Standard Model values of zero. A simultaneous fit to both the 
parameters gives the 95% confidence level limits: 

-0.008 GeV"2 < ao/A^ < 0.005 GeV~^ and - 0.007 GeV"^ < a./A^ < 0.011 GeV"^ , 

as shown in Figure [HI A correlation coefficient of —57% is observed. The experimental sys- 
tematic uncertainties and those on the Standard Model e+e" — >• Z77 qq77 cross section 
predictions are taken into account in the fit. 
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Appendix 

To allow the combination of our results with those of the other LEP experiments, the cross 
sections a are also measured in the more restrictive phase space obtained by modifying the 
conditions (2) and (4) into | cos6'^| < 0.95 and cos6'^q < 0.9, respectively. The results are: 

cT(182.7GeV) = 0.11 ± 0.11 ± 0.01 pb (SM : 0.233 ± 0.003 pb) 

a(188.7GeV) = 0.28 ± 0.07 ± 0.03 pb (SM : 0.214 ± 0.003 pb) 

a(194.5 GeV) = 0.15 ± 0.07 ± 0.02 pb (SM : 0.197 ± 0.003 pb) 

a(200.2 GeV) = 0.15 ± 0.07 ± 0.01 pb (SM : 0.185 ± 0.003 pb). 

The first uncertainty is statistical, the second systematic and the values in parentheses indicate 
the Standard Model predictions. The samples at ^/s = 192 GeV— 196 GeV and ^/s = 200 GeV— 
202 GeV are respectively merged into the -y/s = 194.5 GeV and y/s = 200.2 GeV ones. 
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Figure 1: Diagrams of a) the Standard Model contribution to e"'"e~ — > Z77 signal and 
"non-resonant" background, b) the background from direct radiation of photon from 
the quarks, c) the background from photons, misidentified electrons or unresolved 
TT^s originating from hadrons and d) the anomalous QGC diagram. 
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Figure 2: Distributions of a) the invariant mass Mqq of the hadronic system, b) 
the relativistic velocity of the reconstructed Z boson, c) the energy E^i of the 
most energetic photon and d) the angle uj between the least energetic photon and 
the nearest jet. Data, signal and background Monte Carlo samples are shown for 
= 192 GeV — 202 GeV. The arrows show the position of the final selection 
requirements. In each plot, the selection criteria on the other variables are applied. 
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Figure 4: The cross section of the process e^e~ — Z77 qq77 as a function of the 
centre-of-mass energy. The signal is defined by the phase-space cuts of Equations 
(1) — (4). The width of the band corresponds to the Monte Carlo statistics and theory 
uncertainties. Dashed and dotted lines represent anomalous QGC predictions for 
ao/A^ = 0.015 GeV-2 and Oc/A^ = 0.015 GeV-^, respectively. 
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Figure 5: Distributions for the least energetic photon, a) the energy E^2, b) the 
cosine of its polar angle | cos 9^2 1 , c) its transverse momentum Pt^2 with respect to 
the beam axis. Data, signal and background Monte Carlo are displayed for the full 
data sample together with QGC predictions. 
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Figure 6: Two dimensional contours for the QGC parameters ao/A^ and ac/A"^. 
The fit result is shown together with the Standard Model (SM) predictions. 
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